Stationary states and energy cascades in inelastic gases 
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We find a general class of nontrivial stationary states in inelastic gases where, due to dissipation, 
energy is transfered from large velocity scales to small velocity scales. These steady-states exist for 
arbitrary collision rules and arbitrary dimension. Their signature is a stationary velocity distribution 
f(v) with an algebraic high-energy tail, f(v) ~ v~ a . The exponent a is obtained analytically and it 
varies continuously with the spatial dimension, the homogeneity index characterizing the collision 
rate, and the restitution coefficient. We observe these stationary states in numerical simulations in 
which energy is injected into the system by infrequently boosting particles to high velocities. We 
propose that these states may be realized experimentally in driven granular systems. 
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Energy dissipation has profound consequences in gran- 
ular media. It is responsible for collective phenomena 
such as hydrodynamic instabilities [jj, |2j , shocks |3j, Ij, |5| , 
and clustering 0, 0, II] • The inelastic nature of the col- 
lisions remains crucially important in dilute settings and 
under vigorous forcing where, in contrast with molecu- 
lar gases, there is no energy equipartition 0, Hof and 
the vel ocity distributions are typically non-Maxwellian 
[ill H2I Il3l Hil Hoi Hij ]. In this letter, we show that en- 
ergy dissipation also results in self-sustaining stationary 
states where energy cascades from large velocity scales 
to small velocity scales, and we propose that these states 
may be realized experimentally in driven granular gases. 

Hydrodynamic theory of granular flows is formulated 
using inelastic gases as a starting point [13, [la llil |2Q> 
Elf . Without forcing, dissipation is quantified via the en- 
ergy balance equation, dT/dt = —T, where T — (v 2 ) is 
the temperature and T the dissipation rate. Collisions 
involving a pair of particles with relative velocity Av oc- 
cur with rate cx (Aw) A and result in energy loss cx (Aw) 2 . 
Assuming that a single velocity scale characterizes the 
system, T ~ ((Av) 2+X ) ~ T 1+A / 2 , leads to Haff's cooling 
law, T ~ t~ 2 l x (exponential decay occurs when A = 0). 
Hence, the temperature decays indefinitely and the veloc- 
ity distribution f(v) approaches the trivial steady-state 
f(v) — > S(v) where all particles are at rest, a stationary 
state that can be considered as a dynamical fixed point. 
However, the energy balance equation assumes a finite 
dissipation rate. This need not be the case! 

Our main result is that there is a general class of non- 
trivial stationary states where the velocity distribution 
has an algebraic high-energy tail 
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as v — > 00. First, we obtain this result in one-dimension 
and then, generalize it to arbitrary dimension. 

Our system is an ensemble of identical particles that 
undergo inelastic collisions 



with ui t 2 ("1,2) the pre-collision (post-collision) veloci- 
ties. The collision parameters p and q obey p + q = 1 
so momentum is conserved. In each collision, the rel- 
ative velocity is reduced by the restitution coefficient 
r = 1 — 2p and the energy loss equals pq(u\ — u 2 ) 2 . We 
consider general collision rates \u\ — u 2 \ x with A the ho- 
mogeneity index. For particles interacting via the central 
potential U(r) ~ r~ a then A = 1 - 2^- so A < 1 [22]. 
Hard-spheres, A = JL model ordinary granular media. 
Maxwell molecules 23], A = 0, model granular particles 
with dipole interactions, such as magnetic particles or 
particles immersed in a fluid |24j . 

We seek stationary velocity distributions f(v) that 
obey the Boltzmann equation 



dui du 2 |«i - u 2 \ x f(ui)f(u 2 ) 

x [8(v — pu\ — qu 2 ) — S(v — Mi)] = 0, 



(3) 



reflecting balance of gain and loss due to collisions. 

The case A = is instructive as the velocity dis- 
tribution can be obtained explicitly. The Boltzmann 
equation is a simple convolution because the collision 
rate is constant, and it is conveniently studied using 
the Fourier transform F(k) = J dvf(v)e lkv that satis- 
fies the nonlinear, nonlocal equation F(k) — F(pk)F(qk). 
Normalization implies F(0) = 1. For elastic collisions 
{p = 0) any distribution is stationary, but this is just a 
one-dimensional pathology. Generally, for any value of 
the restitution coefficient, there is a family of solutions 
F(k) = exp ( — t>o|&|) with arbitrary typical velocity Vo- 
As a result, the velocity distribution is Lorcntzian 
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We note that this distribution is independent of the resti- 
tution coefficient. The typical velocity vq is finite, but the 
average energy is infinite. 

This distribution does not evolve under the inelastic 
collision dynamics and in particular, the energy density 
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remains constant. From the energy balance equation, one 
might have expected that the velocity distribution con- 
stantly narrows because particles dissipate energy, but 
quite surprisingly, the heavy tail of the velocity distribu- 
tion acts as a heat bath, balancing the energy dissipa- 
tion and maintaining a steady velocity distribution. We 
conclude that, in addition to the trivial state where all 
particles are at rest, there is another fixed point, a self- 
sustaining stationary state. 

For general collision rules, the characteristic exponent 
a can be found analytically. For large v, the convolution 
in Eq. (JSJ is governed by the product f(u\)f(u2) with one 
of the pre-collision velocities large and the other small be- 
cause the distribution decays sharply at large velocities. 
The Boltzmann equation includes a loss term and a gain 
term. For the loss term, u\ — v is large and ui is small. 
For the gain term, there are two separate contributions. 
Either v = pu\ and then ui is small or v — qui and then 
Mi is small. In either case, the double integral separates. 
The integral over the smaller velocity equals one. With 
the remaining integral over the larger velocity, denoted 
by u, the nonlinear Boltzmann equation becomes linear 
for large v 



du\u\ x f(u) [S(v-pu)+S(v-qu)-5(v-u)] = 0. (5) 



Consequently, the tail of the distribution satisfies the 
functional equation 
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describing cascade of energy from large velocities to small 
ones, v — > (pv,qv). The power-law decay Q satis- 
fies this equation when p <T ~ 1 ~ x + q cr ~ 1 ~ A = 1, and since 
p + q = 1, the characteristic exponent is 



a = 2 + A. 



(7) 



In one-dimension, a is independent of the restitution co- 
efficient. 

We comment that the full nonlinear Boltzmann equa- 
tion assumes molecular chaos as the two-particle veloc- 
ity distribution is a product of one-particle distributions. 
However, the equation governing the tail is linear and 
it is valid under less stringent conditions. The only re- 
quirement is that energetic particles are weakly corre- 
lated with slower ones. 

In arbitrary dimension d, the collision law is 



vi = ui - (1 -.p)(ui - u 2 )-nn. 



(8) 



Momentum transfer occurs only along the normal direc- 
tion with n a unit vector parallel to the impact direc- 
tion connecting the particles. The energy dissipation 
equals p(l — p)|(ui — U2)-n| 2 and the collision rate is 
(ui -u 2 )-n| A . 



We derive the linear equation governing the tail 
of the distribution directly from the collision rule. 
With the large velocity u, its magnitude u = |u| and 
/i=(u-n) 2 , the collision rate is u A /W 2 . The cas- 
cade process is v^>(av,(3v) with the stretching pa- 
rameters a = (1 — p)/! 1 / 2 and (3 = [1 — (1 — p 2 )^] 1 / 2 for 
v = (1 — p)u • fin and v = u — (1 — p)u ■ fin, respec- 
tively [23. Integrating over the impact angle n and the 
velocity magnitude u, Eq. (0) generalizes as follows 



du 



i A+d - 1 /(") [S(v-au)+S(v-pu) - /(w)]/i A/2 ^ = 0. 



Here, the brackets are used as shorthand for the angular 
integration, (g) = J dhg(h). Integrating over u yields 
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A / 2 ) = 0. (9) 



The power-law decay satisfies this equation pro- 
vided that the exponent a is root of the equa- 
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gular integration is performed using spherical coordi- 
nates. Given /x = cos 2 6, then dn oc (sm8) d ~ 2 d9, and 
(g) oc f Q d/i/i~2 (l — The equality for the ex- 

ponent er can be written in terms of the gamma function 
and the hypergeometric function [2fij | 
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In the cascade process v — » (av,0v), the total veloc- 
ity magnitude increases (a + (3 > 1) even though the to- 
tal energy decreases (a 2 + f3 2 < 1). Utilizing these in- 
equalities, we deduce the bounds l<er — d — X < 2 from 
the embedded equation for the characteristic exponent. 
The lower bound (J7J) is realized in one-dimension where 
collisions are head-on (/j, = 1). The upper bound is ap- 
proached, a — > d + 2 + A, in the quasi-elastic limit, r — > 
1. This, combined with the known Maxwellian distribu- 
tion that occurs for elastic collisions, demonstrates the 
singular nature of the quasi-elastic limit. Interestingly 
the energy may be either finite or infinite, depending on 
whether a is larger or smaller than d + 2. In either case, 
the integral underlying the dissipation rate is divergent, 
a general characteristic of the stationary distributions. 

Generally, the exponent a increases monotonically 
with the spatial dimension d, the homogeneity index 
A, and the restitution coefficient r. For completely in- 
elastic (r = 0) hard spheres, we find a = 4.14922 and 
a = 5.23365 for d — 2, 3, respectively. The exponents for 
completely inelastic Maxwell molecules are a — 3.19520 
and a = 4.28807 for d = 2, 3. These values provide lower 
bounds on a with respect to r, as shown in figure 1. 

A remarkable feature is that the characteristic expo- 
nent changes continuously with the parameters d, A, and 
r. Hence, the tails are not universal. The power-law 



FIG. 1: The exponent a as a function of the restitution co- 
efficient r for A = (circles) and A = 1 (squares). The lower 
curves correspond to d — 2 and the upper to d = 3. 
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FIG. 2: The velocity distribution f(v) versus v for Maxwell 
molecules in one-dimension. The velocity was rescaled to 
compare with an ordinary Lorentzian. 

decay stands in contrast with the stretched exponen- 
tials, f(v) ~ exp(— v v ), with v = A and v = 1 + A/2, 
respectively, for unforced and thermally forced inelastic 
gases 0, |3 EH E3- Previously, algebraic tails with 
different exponents were found but for non-stationary 
distributions describing freely cooling Maxwell molecules 

null. 

These steady-state distributions can be realized, up to 
some cutoff, in finite systems of driven inelastic parti- 
cles. The key is to inject energy at a very large velocity 
scale. Numerically, we used the following simulation. Ini- 
tially, we start with an innocuous velocity distribution, 
uniform with support in [—1 : 1]. Particles collide in- 
elastically, and an energy loss "counter" keeps track of 
the aggregate energy loss. With rate 7, small compared 
with the collision rate, a particle is selected at random 
and it is energized by an amount equal to the aggregate 




FIG. 3: The velocity distribution f(v) versus v for Maxwell 
molecules (a) and hard spheres (b) in one-dimension (upper 
lines) and two-dimension (lower lines). The theoretical pre- 
diction is shown for reference. 

energy loss. Subsequently, the energy loss counter is re- 
set to zero. The rationale behind this simulation is that 
the total energy remains practically constant and more 
importantly, that energy is injected only at the tail of the 
distribution. In effect, this procedure does not alter the 
purely dissipative dynamics, except for setting a scale for 
the most energetic particles. 

We simulated completely inelastic hard-spheres and 
Maxwell molecules in one- and two-dimensions. We used 
N = 10 7 (N = 10 5 ) particles and injection rate 7 = 10~ 4 
(7 = 1(T 2 ) for A = (A = 1). We verified that: (i) 
after a transient, the velocity distribution becomes sta- 
tionary, (ii) the velocity distribution is Lorentzian for 
one-dimensional Maxwell molecules (figure 2), (iii) the 
velocity distribution has an algebraic tail, and (iv) a is 
consistent with theory (figure 3). 

The very same stationary distributions can be reached 
using other simulation methods. For example, particles 
may be re-energized such that their velocity is drawn 
from a Maxwellian distribution with a typical energy pro- 
portional to the system size (the data presented for one- 
dimensional hard-spheres are from such a simulation). 
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We conclude that the simulations confirm the existence 
of the nontrivial stationary states with power-law tails. 
These steady-states are stable fixed points as the system 
is driven into them even when starting from compact 
distributions. Moreover, stability analysis of the time 
dependent version of Eq. © shows that the stationary 
distribution (JIJ is stable with respect to perturbations 
consisting of steeper algebraic tails. 

Clearly, if f(v) is a steady-state so is VQ d f(v/vo) for 
an arbitrary typical velocity Vq. The injection protocol 
selects wo- Suppose that particles are boosted at rate 
7 per particle to velocity V, a scale that sets an upper 
cutoff on the velocity distribution j^Oj. The energy in- 
jection rate is jV 2 , and the energy dissipation rate is 
T ~ (v 2+x ) with (g) = f V dvv d ~ 1 f(v)g(v). Energy bal- 
ance, r ~ jV 2 , relates the injection rate, the injection 
scale, and the typical velocity scale, 7 ~ V x {V/vo) d ~ a . 

In our simulations, energy is maintained constant, 
(v 2 ) ~ 1. When a < d+2, the constant energy constraint 
implies y d + 2 - a ^ VQ~ a , that combined with energy bal- 
ance reveals how the maximal velocity, V ~ 7 _1 /(2-A) j 
and the typical velocity, vq ~ 7 <"- d )( 2 -*) , scale with the 
injection rate. When a > d + 2, the initial conditions set 
the typical velocity: (v 2 ) ~ v 2 implies vq ~ 1, and energy 
balance yields V ~ 'y-VC 0- — The simulations are 
consistent with these scaling laws. For example, V « 10 2 
and v ps 10~ 2 for one-dimensional Maxwell-molecule 
simulations with 7 = 1CP 4 . This, combined with the 
simulations, demonstrates that finite energy and infinite 
energy cases differ quantitatively but not qualitatively. 

In conclusion, we have shown that inelastic gases have 
nontrivial stationary states where the velocity distribu- 
tion has an algebraic high-energy tail. These steady- 
states are generic in that they exist for arbitrary colli- 
sion rates, arbitrary collision rules, and arbitrary spatial 
dimensions. Typically, energy is injected at large veloc- 
ity scales, and it is dissipated over a range of smaller 
velocities. Qualitatively, this energy cascade picture is 
analogous to fluid turbulence. However, without energy 
conservation, there is an important difference with the 
conventional Kolmogorov spectra: the characteristic ex- 
ponents are not universal and they do not follow from 
dimensional analysis |3lj . 

In an infinite system, there is perfect balance between 
collisional loss and gain and the high-energy tail enables 
the stationary state, while in finite systems, energy in- 
jection maintains the steady-state. The latter is relevant 
to granular gas experiments that are typically performed 
in steady-state conditions. We propose that stable dis- 
tributions and energy cascades may be realized using an 
experimental set-up in which energy is added at large 
velocity scales via rare but powerful energy injections. 
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